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Abstract The compatibility of the fast tachocline scenario with a flux transport 
dynamo model is explored. We employ a flux transport dynamo model coupled 
with simple feedback formulae relating the thickness of the tachocline to the 
amplitude of the magnetic field or to the Maxwell stress. The dynamo model 
is found to be robust against the nonlincarity introduced by this simplified fast 
tachocline mechanism. Solar-like butterfly diagrams are found to persist and, 
even without any parameter fitting, the overall thickness of the tachocline is 
well within the range admitted by hclioseismic constraints. In the most realistic 
case of a time and latitude dependent tachocline thickness linked to the value 
of the Maxwell stress, both the thickness and its latitude dependence are in 
excellent agreement with seismic results. In the nonparamctric models, cycle 
related temporal variations in tachocline thickness are somewhat larger than 
admitted by helioseismic constraints; we find, however, that introducing a further 
parameter into our feedback formula readily allows further fine tuning of the 
thickness variations. 

Keywords: dynamo, tachocline 
1. Introduction 

Flux transport dynamos are the most widely discussed scenario for the origin of 
solar activity. In these models, the a-effect responsible for toroidal — » poloidal 
flux conversion is identified with the tilting of the axis of active regions relative 
to the azimuthal direction. The poloidal field generated by this a-effect concen- 
trated near the surface is then transported by the meridional circulation to high 
latitudes and then down to the bottom of the convective zone (Babcock-Lcighton 
mechanism) . The return branch of a one-celled circulation flow pattern will then 
advect the field back towards the equator while differential rotation amplifies 
it to a toroidal field of increasing strength (Q-effect). When the toroidal field 
strength reaches a critical limit, flux emergence driven by the Parker instability 
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will give rise to active regions on the surface, with an axis tilted away from 
the azimuthal plane due to the Coriolis force, thereby closing the cycle (see 
e.g. Charbonneau, 2010 for a more detailed discussion). An attractive feature of 
this dynamo model is that it not only models the regular features of the solar 
cycle (Choudhuri, Schiissler, and Dikpati, 1995; Dikpati and Charbonneau, 1999; 
Chatterjce, Nandy, and Choudhuri, 2004) but also reproduces many irregular 
features of solar cycle (Charbonneau and Dikpati, 2000; Choudhuri and Karak, 
2009; Karak, 2010; Karak and Choudhuri, 2011; Choudhuri and Karak, 2012). 

In order to avoid excessive buoyant loss of toroidal flux, the fi-effect must 
take place in a stably stratified layer just below the convective zone. Helio- 
seismic inversions (Charbonneau et al., 1999; Basu and Antia, 2001) indicate 
that this layer coincides with the solar tachocline, i.e. the thin transitional layer 
between the rigidly rotating interior and the differentially rotating convective 
zone. The narrow radial extent of this layer (scale depth of 5-10 Mm or less) 
implies a strongly anisotropic angular momentum transport mechanism, much 
more effective in the horizontal direction than in the radial direction. Spiegel 
and Zahn (1992) have shown that this is possible by taking strong anisotropic 
turbulent viscosity in the horizontal direction. However several authors (e.g., 
Gough and Mclntyre, 1998; Dikpati and Gilman, 1999) have found caveats in 
this purely hydrodynamical model and improved by including additional physics. 
Other most plausible candidate for such anisotropic momentum transfer is the 
Maxwell stress in a predominantly horizontal magnetic field configuration. One 
school of thought in this stream is that a weak fossil magnetic field in the 
radiative zone may solve this problem (Rudigcr and Kitchatinov, 1997; Gough 
and Mclntyre, 1998; MacGregor and Charbonneau, 1999). Unfortunately the 
simulations (Garaud, 2002; Brun and Zahn, 2006; Strugarek, Brun, and Zahn, 
2011) of this so-called slow tachocline model are not able to provide uniform 
rotation in the radiative zone even by including penetrating flows from the 
convection zone, such as plumes or meridional circulation to confine the fossil 
magnetic field (however Garaud and Garaud, 2008 succeeded in this problem). 
Other school of thought is that as the dynamo generated toroidal field resides 
in the tachocline, it is plausible to assume that this field is responsible for the 
confinement of the tachocline. The feasibility of this so-called fast tachocline 
scenario has been demonstrated by Forgacs-Dajka and Petrovay (2001, 2002) 
and Forgacs-Dajka (2003). 

In the fast tachocline scenario, the thickness of the tachocline depends on 
the magnetic field. This thickness, on the other hand, is one input parameter 
of flux transport dynamo models and thus it may be expected to influence the 
amplitude and configuration of the magnetic field. This coupling introduces a 
nonlincarity into the dynamo that may affect both the dynamo and the structure 
in a number of ways. Indeed, it is not a priori clear whether the fast tachocline 
scenario and flux transport dynamos are compatible at all, i.e. whether they 
can give rise to a finite amplitude oscillatory field with a tachocline of realistic 
thickness and with cyclic variations within the observational bounds. 

The aim of the present paper is to investigate this issue. In a completely self- 
consistent approach one should couple the angular momentum transfer equation 
to the induction equation solved in dynamo models. Instead of this more generic 
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approach, as a first step, here we parameterize the dependence of tachocline 
width d t on the poloidal field amplitude B p with simple algebraic formulae that 
still retain the essential physics of the problem. Details of our problem setup are 
presented in Section [2] Section [3] presents the results which are then confronted 
with observations in Section 4. Section [5] concludes the paper. 



2. Problem setup 

2.1. Flux transport dynamo model 

An axisymmetric magnetic field can be represented in the form 

B = B(r, 9)e ( j } + V x [A(r, 9)e <jl ], (1) 

where B(r,9) and A(r,9) respectively correspond to the toroidal and poloidal 
components. Then the evolution of magnetic fields in the flux transport dynamo 
model arc governed by the following two equations: 
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with s = r sm9. 

Here v is the meridional flow which has the following analytical form. 
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factor which determines the maximum value of the latitudinal component of the 
meridional circulation v$. 

77 is the turbulent magnetic diffusivity which has the following form: 



V( r ) = VRZ 



Vscz 



1 + erf 2 



,r -r B cz \ 
dt J 



1 + erf 
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(6) 



with r BC z = 0.7R, d t = 0.03R, d 2 = OMR, r surf = 0.95i?, r] RZ = 5 x 10 s 
cm 2 s _1 , 77SCZ = 5 x 10 10 cm 2 s _1 and 7j sur f = 2x 10 12 cm 2 s" 1 . 
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Figure 1. Reference model with fixed tachocline thickness dt- Contours show the butterfly 
diagram of the toroidal field in the tachocline, the field strength being 0.91 X 10 5 G at the 
innermost contours and half that value at the outermost contours. Blue contours correspond to 
positive toroidal field whereas red contours correspond to negative toroidal field. The greyscalc 
background shows the weak diffuse radial field on the solar surface. 



S(r, 9; B) is the coefficient which describes the generation of poloidal field at 
the solar surface from the decay of tilted bipolar sunspots (Babcock-Lcighton 
process). It has the following form 



S(r,0;B) = 



a{r, 9) 



l + (B(r t ,9)/B Q y 



■B{r u 



(7) 



where a(r, 0) = 2-2. 
r 4 = Q.95R, r 5 = 



1 + erf 



1 - erf 



R, d 4 = 0.05R, d 5 = OMR, 7 
B = 4 x 10 4 G. 

f2 is the solar rotation which has the following form 
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where r t (the position of the tachocline) = 0.7-R, dt (the thickness of the tachocline) 
= 0.03i?, H R z/2tt = 432.8 nHz, Cl s cz(9) = Oeq + a 2 cos 2 {9) + a 4 cos 4 (9), with 
E q/27t = 460.7 nHz, a 2 /27r = -62.69 nHz and a 4 /2ir = -67.13 nHz. 

With the appropriate boundary conditions we solve Equations (2) and (3) in 
a full sphere of the meridional plane with 0.55i? < r < R, < 9 < it to study 
the evolution of magnetic fields (see Chatterjee, Nandy, and Choudhuri, 2004 
for details). 

Once the transient state is over the butterfly diagram looks like the one shown 
in Figure [TJ In this figure the contours show the toroidal field in the tachocline 
and the shaded background shows the strength of the radial field on the solar 
surface. 

In the next section, we shall discuss the dependence of the tachocline thickness 
on the magnetic field based on the fast tachocline model and then we shall study 
its effect on this flux transport dynamo model. 
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2.2. Tachocline depth parametrization 

Tachocline thickness is introduced in various ways by various authors, so care 
must be taken when comparing results from different papers. For a comparison 
it is useful to define the relative residual angular velocity £1 as 

"SCZlyj — "RZ 

Clearly, Cl takes the value in the radiative zone, f at the top of the tachocline. 
The radius r t of the tachocline is defined by 0(rt, 9) = 0.5. 

If O varies by a factor F in the radius range \r t — x/2, r t + x/2], then clearly 
x/H = \nF where H is the scale height of the tachocline. From Equation ([5]) 
it then follows that our tachocline thickness parameter dt = 1.65H. Similarly, 
for the tachocline thickness parameter w used by Antia and Basu (2011) in 
their helioseismic study of tachocline properties we find w = 0.39H, i.e. dt — 
4.3iu. This relation must be taken into account when comparing our results to 
observational constraints. 

An approximate analytic formula relating the mean tachocline scale height 
H to the amplitude of an imposed poloidal field in a periodic shear layer was 
derived by Forgacs-Dajka and Petrovay (2001). It reads 

2 Prrirju (1 + r,/uH 2 )(l + Pr V /ujH 2 ) 
p H 2 1 + 2Yvi 1 /loH 2 ' 1 J 

where uj = 2n/22 years is the dynamo frequency, Pr is the Prandtl number, 
i] is the turbulent magnetic diffusivity, r t is the radius of the tachocline and 
V p is the amplitude of the oscillatory poloidal magnetic field in Alfvenic units. 
In a turbulent medium the effective Prandtl number is expected to be order 
of unity, while its exact value is unknown. Therefore, in what follows we shall 
simply substitute dt for H in Equation (jlOp . ignoring the factor 1.65 between 
these scales as an appropriate choice of Pr can always offset this factor anyway. 

Figurc[5]presents the variation of B p (= V p (47rp) 1 / 2 with p = 0.1 g/cm 3 ) based 
on the relation (| 10[) for three different values of the diffusivity. Green curves are 
simpler analytical fits to each curve with a function of the form B p = Crj/d 2 or 

= % (11) 



where C = 6 x 10 10 G-s, and B p = \J {B 2 + Bf). 

In a first attempt to relate tachocline thickness to the dynamo field by a simple 
power law formula, in Section 13.11 we shall employ formula (jllj) as a guidance 
in our choice of the exponent. It should be emphasized that, as Equation (jTTJ) 
refers to cycle averaged quantities, it should not be viewed as rigorous derivation 
of our coupling formulae (|2"Tj) and ([231) ). Rather, it provides a distant analogy 
guiding us in our choice of a physically motivated exponent value. In particular, 
in the fast tachocline scenario the tachocline is confined by the Maxwell stress 
B p B which also involves the toroidal field strength B. The reason why B docs 
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Figure 2. Relation between tachocline width dt and poloidal field amplitude B v for three 
different values of magnetic diffusivity, r] = 10 12 (solid), tj = 3 X 10 11 (dotted), r] = 10 11 
(dashed). Green curves show the function B p = 6 X 10 10 r;/d 2 . 

not explicitly appear in Equation is that the toroidal held is generated by 
the windup of the poloidal field, so its overall amplitude is ultimately determined 
by Bp. However, this is clearly only true for the overall amplitudes and not for 
the actual values of B p and B at any given point and instant: so the simple 
approximate relationship is, strictly speaking, only valid for the (latitudinal 
and temporal) average field amplitude and average tachocline thickness and its 
application to relate dt and B p at any given point in space and time is not well 
justified. 

For this reason, in Section 13.21 we shall try to further improve on the model 
by relating dt directly to the Maxwell stress. We recall from Forgacs-Dajka and 
Petrovay (2001) that the simple model giving rise to Equation (jTTJ) consists 
of a semi- infinite layer bounded from above at depth z = ("bottom of the 
convective zone" ) where a periodic shear flow is imposed in the y direction: 

VyO = VQCOs(kx). (12) 

An oscillatory horizontal "poloidal" field is prescribed in the x direction through- 
out the volume as 

B x = B p cos(ut). (13) 
Introducing v = v y and using Alfven speed units for the magnetic field 

V p = BpiA-Kp)- 1 ' 2 b = ByiAnp)- 1 / 2 , (14) 
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the azimuthal component of the equation of motion reads 

dtv = V p cos(u;t)d x b — vV 2 v, (15) 

where v denotes viscosity. Solutions may be sought in the form 

v = v(x,z) + v'(x,z)f(ui), b = b'(x,z)f(ojt + 4>), (16) 

where / is a 27r-periodic function of zero mean and of amplitude 0(1). (a denotes 
time average of a, while a' = a — a.) The (temporal) average of Equation (I15|) 
then reads 



= V p cos(cjt)f(u)t + 4>)d x b' - vV 2 v. (17) 
For an order of magnitude estimate of the terms of this equation, we assume 



cos(uit) f (uit + 4>) 6 0(1) (i.e. no "conspiracy" between the phases, a rather 
natural assumption). As H <C R we may approximate V 2 ~ H~ 2 . Estimating 
the other derivatives as dt ~ w and d x ~ -R -1 , (fT7|> yields 

V p b'/R~ isv/H 2 , (18) 



vRv 

1 ^ 



H - w (19) 



or, switching back to conventional magnetic field units and from H to d- 



« - 0" < 20 ' 

Plugging in the correct dimensional values we find C ~ 2 x 10 15 G 2 -s — this 
is the value we use in Section 13.21 below where, again, we shall employ formula 
(|20[) as a guidance in our choice of the exponent of our simple feedback formula. 
While the link between tachocline thickness and Maxwell stress is more direct 
than between d t and B p alone, it should still be kept in mind that equation 
(1201) refers to cycle averaged quantites, so it should not be viewed as a rigorous 
derivation of our coupling formulae (|26l) and (|29[) . Rather, it provides an analogy 
guiding us in our choice of a physically motivated exponent value. 



3. Results 

3.1. Tachocline thickness linked to poloidal field amplitude 

In the fast tachocline scenario the thickness of the tachocline may vary as a func- 
tion of both time and latitude. For a fully consistent treatment of this variation 
the equation of motion should be coupled to the dynamo equations. However, in 
the present preliminary exploration of the problem we are only interested in the 
general stability properties and robustness of flux transport dynamo against the 
kind of nonlinearity introduced by the fast tachocline mechanism. Therefore, we 
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simplify the problem by an appropriate parametrization of the dependence of dt 
on the magnetic field. For this purpose any arbitrary numerical relation between 
tachoclinc thickness and magnetic field may do as long as it results in dt values 
of the right order of magnitude and it has the expected property of resulting in 
a thinner tachocline for stronger magnetic fields. 

The numerical relation we first consider here is Equation (HID . As we already 
remarked in Section 2.2, by its derivation, the formula (|11[) given above only 
relates the cycle and latitude averaged mean value of the tachocline thickness 
to the amplitude of the variation of the poloidal field, and our use of it to link 
local and momentary values of these variables is arbitrary. Nevertheless, by its 
construction, Equation (|11D does have the required properties of resulting in d t 
values of the right order of magnitude and in a thinner tachocline for stronger 
magnetic fields, so it will suffice for the purpose of a first exploration. 

For a first study we shall not consider any latitudinal variation of the tachoclinc. 
Substituting the appropriate turbulent value for 77 in the tachocline we then have 

di = (21) 
' B p {t) V ' 



where 



1 



r t +d 



m = ^ / V(r) dr (22) 



rt—ao 



with d = 0.015i?. 

9 /-t/2 



B p (t) = - [ B p {9,t)d6 (23) 
n Jo 

is the latitude averaged poloidal field amplitude while B p (0, t) is the local radial 
mean value of the poloidal field calculated as 

1 rrt+do 

B P (0,t) = — B p {r,6,t)dr. (24) 

M J r ._ dn 



(Recall that B p = y/ (£>,,? + Bg), i.e. the modulus of the poloidal field. This 
ensures that B p , as given in the above equations, remains positive at all times.) 

Now, starting from the relaxed state of our (fixed dt) model discussed in 
Section 2.1, we allow dt to vary at every time step of our simulation, calculating 
B p (t) as above to get the value of dt from Equation (j2Tj) . The calculation is run 
again for several solar cycles until it relaxes to a nearly steady cyclic behavior. 
In Figure 02 we show a few solar cycles long clip just after this modification. 
It is noteworthy that this kind of tachoclinc structure calculated from the fast 
tachocline model does produce a periodic solar-like solution. In the upper panel, 
the butterfly diagram of toroidal and radial field components is shown. The 
appearance of this plot is quite similar to the one found in Figure [1] in Sec- 
tion 2.1. Therefore we may conclude that the fast tachoclinc model and the 
flux transport dynamo model are compatible and a nonlinear coupling between 
tachocline thickness and field amplitude does not lead to a breakdown of the 
solar-like solution. 
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Figure 3. Upper panel: same as Figurc^but for variable tachocline thickness dt as given by 
cq. (1 2 1 1) . Lower panel: variation of tachocline thickness dt with time. 

The lower panel shows the variation of dt with solar cycle. The amplitude of 
this variation is about ±25 %. 

Next, we try to explore the latitudinal structure of the tachocline by admitting 
d t to depend on 6, i.e. instead of Equation ([2~Tj) we use 



d 2 t = "=~77T~T- (25) 

Figure |4] presents the result for this calculation. Upper panel shows the butter- 
fly diagram which still resembles most characteristics of the observed butterfly 
diagram. 



3.2. Tachocline thickness linked to the Maxwell stress 



In our next parametric study we link the tachocline thickness dt to the Maxwell 
stress by formula (|20|) . Again, for our first study we disregard the latitude 
dependence in d t : 

d l = _ C '^_ (26 ) 
4 B p (t)B(t) V ' 
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Figure 4. Upper panel: same as Figure^ but for variable tachoclinc thickness dt as given 
by eq. il 2 -"ill . Lower panel: variation of tachoclinc thickness dt with time at different latitudes. 
Dash-dotted, solid, dashed and dotted lines are the values at 75°, 60°, 45° and 15° latitudes, 
respectively. 



The mean values of B are here defined in a manner analogous to those of B p : 

B(t) = - / B(0,t)d6 (27) 
n Jo 

is the latitude averaged toroidal field amplitude while B(9,t) is the local radial 
mean value of the toroidal field calculated as 



^ e ^ = ^T ° \B{r,6,t)\dr. 



(28) 



(Again, using the modulus of the toroidal field in the integrand ensures that B, 
as given in the above equations, remains positive at all times.) 
The result is shown in Figure [5] 

Finally, to explore the latitude dependence of tachocline thickness we set 



C% 



B p {6,t)B(6,t) 



(29) 



This yields what we may regard our most realistic model. Figure |6] shows the 
result. 
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Figure 5. Same as Figure[3]but here the tachocline thickness dt is determined by eq. (I26II . 
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Figure 6. Same as Figure [4] but here the tachocline thickness dt is determined by eq. II29H 
(most realistic model). In the lower panel, the dash-dotted, solid, dashed and dotted lines are 
the values of dt at 75°, 60°, 45° and 15° latitudes, respectively. 
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Figure 7. Variation of tachocline thickness dt (a); of the amplitude variation of the tachocline 
thickness A/dt (b); and of toroidal field strength at 15° latitude (c) as functions of ft (the 
exponent in Equation 1301 



4. Discussion 

The properties of the variable tachocline resulting from our models may be 
compared to empirical (hclioscismic) constraints. An in-depth study of tachocline 
thickness as a function of latitude and time was recently performed by Antia 
and Basu (2011) (also see Basu and Antia, 2001). Converting their tachocline 
thickness parameter w to our parameter dt as discussed in Section 2.2 above, 
their finding was that d t / R varies from approximately 0.02 to 0.1 with increasing 
hcliographic latitude. Temporal variations of up to ±50 % are observed; however 
these are dominated by variation on the time scale of a few years. (Variations 
on such time scales may be produced by shear instabilities, cf. Micsch, 2007) On 
the longer time scale of 11 years no systematic variations were detected; from 
visual inspection of the relevant figures, this sets an upper limit of ~ ±20 % on 
cycle related changes in tachocline thickness. 

The mean value of tachocline thickness in Figures [3H6J lies around 0.03i?. In 
the most realistic case, Figure [51 the thickness varies from 0.02 to 0.1 as we 
move from low to high latitudes, in perfect agreement with the observations. 
This demonstrates that the simple analytical feedback formula (|29p not only 
gives a solar like dynamo solution in the flux transport dynamo model but 
also reproduces the realistic value and the observed latitude dependence of the 
tachocline thickness. 

The cycle variation of tachocline thickness found in our models is quite sig- 
nificant, up to ±30%. This is somewhat higher than what the observational 
constraints suggest. 

In order to explore how sensitive our quantitative results are to details of the 
feedback formula, we generalize Equation ([2^)1 as 

_ (c%y^ 

dt - [B p (o,t)B{o,t)r m 

Clearly, the 0.5 returns Equation ([2"0|). While using other values of re 

has no clear physical justification, it offers a way to explore the sensitivity of our 
quantitative results to details of the feedback formula chosen. 

We repeat the same calculation at different values of re ranging from 2 to 8. 
Figure [7] presents the results. Panel (a) shows the variation of dt as a function 
of re. Note that as we vary re, we also vary C" in order to fix the mean value of 
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dt at around 0.03i? which is fairly close to the observed value. (As C involves 
the unknown quantity Pr , its value is not well constrained anyway.) Figure [7Jb) 
shows the amplitude A of the departure of d t from its mean value as a function 
of k. A smooth increase with k is found, so lower k values result in rather 
more subdued thickness variations that may be compatible with more stringent 
empirical limits. Following Choudhuri (2003) we should keep in mind that as 
we are using mean field dynamo equations, the Maxwell's stress B p B is consist 
of the mean values of B p and B and thereby these are non-zero only inside the 
flux tubes. Now if (B p )it and (-B)ft are the values of these quantities inside flux 
tubes and / is the filling factor, then we have B p = f(B p )it and B = f{B)it (we 
assumed the same filling factor for both components for the sake of simplicity). 
Therefore it is easy to show that the mean stress appeared in all the calculations 
of Section 3.2 should be replaced by the f(B p )it (-B)ft = B p B/ f However this 
factor / can be absorbed in the C because C is not a strict constant. 

5. Conclusion 

We have explored the compatibility of the fast tachocline model with the flux 
transport dynamo model. For this purpose we employed two simple feedback 
formula, Equations (j25|) and ([29)) , relating the thickness of the tachocline to the 
poloidal magnetic field strength. While the second of these formulae is physically 
more sound, both reflect the expected form of the relationship between cycle 
averaged tachocline thickness and the amplitude of the magnetic field, and they 
do have the required properties of resulting in dt values of the right order of 
magnitude and in a thinner tachocline for stronger magnetic fields. Both cases 
were studied with and without latitude dependence in dt- 

The flux transport dynamo model we considered proved to be robust against 
the nonlinearity introduced by this simplified fast tachocline mechanism. Solar- 
like butterfly diagrams were found to persist and the overall thickness of the 
tachocline is well within the range admitted by helioscismic constraints even 
without any parameter fitting. In the most realistic case of a time and latitude 
dependent tachoclinc thickness linked to the value of the Maxwell stress, d t and 
its latitude dependence are in excellent agreement with seismic results. 

One characteristic of fast dynamo mechanism is a marked cycle dependence 
in tachocline properties. Indeed, without parameter fitting, in all of our models 
the cycle related changes in tachocline thickness are somewhat larger than the 
maximal such variation admitted by helioseismic constraints. On the other hand, 
an exploration of the parameter space indicates that in a slightly generalized form 
of our feedback formula it is possible to find parameter combinations where the 
cycle variation remains within the observational bounds while in other respects 
the good agreement of our model with observations is preserved. 

While there is no straightforward physical justification for such a generalized 
feedback formula, our parametric study shows that results like the compatibility 
of the fast dynamo and flux transport dynamo mechanism, the overall thickness 
of the tachocline or the amplitude of the toroidal field are more robust than 
details like the time and latitude dependence. Physical effects not considered 
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in our model may, then, potentially be held responsible for the deviation of 
the feedback formula exponent from the value of 0.5 derived from a simplified 
analytic model. 

Such further effects may include tachocline instabilities, a change in the 
geometric setup of our dynamo model or a combination of the two. One clear lim- 
itation of the model is that differential rotation is quite strong at high latitudes, 
so one may expect that this is the region where strong toroidal held is produced. 
However, sunspots appear at low latitudes only. In the Surya family of models 
on which our model is based, this problem is bypassed by considering a deeply 
penetrating meridional flow (Nandy and Choudhuri, 2002). While this may be 
unphysical, it has the desired effect of allowing the toroidal field generated at 
higher latitudes to be stored and amplified further until it reaches the lower 
latitudes where sunspot eruptions happen. A promising, and more physically 
consistent way out of this conundrum was suggested by Parfrey and Menou 
(2007) who showed that at latitudes higher than 37° the magnetorotational 
instability (MM) is present in the tachocline, presumably resulting in stronger 
turbulence, a thicker tachocline and a less organized field structure. All this 
just serves to illustrate that, independently of the choice of our feedback for- 
mula, there may be more physical effects to consider before we can hope to 
simultaneously reproduce the correct temporal and latitudinal structure of the 
tachocline. 

A further obvious limitation is that only one particular dynamo model was 
considered. It is well known that properties of flux transport dynamos are quite 
sensitive to variations in the input parameters and assumptions (e.g. Dikpati 
et al., 2002; Guerrero and de Gouvcia Dal Pino, 2007; Jiang, Chatterjee, and 
Choudhuri, 2007; Yeates, Nandy, and Mackay, 2008; Karak and Choudhuri, 2012; 
Karak and Nandy, 2012). It is therefore desirable to extend such studies to other 
dynamo configurations. 

It is nevertheless clear that the use of simple feedback formulae of any form via 
a simplified approach to the problem can only be considered as a first exploration. 
In any case, the completely consistent treatment of the problem should involve 
the full solution of the equation of motion, coupled with the dynamo equations. 
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